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A large nrnnber of Mars craters exist that are nearly tangential to other 
craters. They occur in clusters or as isolated crater doublets. Results of 
probability calculations and a Monte Carlo cratering simulation model show con- 
clusively that many of the Mars craters could not have resulted from random 
single body impact. However, clusters and doublets could be caused by meteoroid 
breakup resulting from stresses induced in the meteoroid by the gravitational 
field of Mars. Calculations are provided that show the required relationships 
between mass, entry angle, impact velocity, and tensile strength of the meteoroid, 
and the separation distance between the centers of the resulting craters relative 
to crater sizes. It is concluded that under certain conditions, doublets should be 
produced on Mars as a direct result of breakup of an impacting meteoroid. The 
impact process can yield nonrandom crater distributions that should be observed 
in different degrees of development on different planetary surfaces. Thus, dis- 
tributions of large craters on different planetary bodies cannot be compared 
directly for interpretations of the geologic history of one planetary surface 


relative to another. 



INTRODUCTION 


Examination of Mariner 6 and 7 photographs has revealed a large number of 
craters on Mars that are tangential or nearly tangential to one another. The 
hypothesis that all craters were formed by random single-body impact was tested 
through examination of the results of a series of Monte Carlo computer simula- 
tions of random impact cratering and through an independent probability compu- 
tation. Both the Monte Carlo simulation and the probability calculation show that it is 
probable that most of the Mars doublets have not been produced by a random process . 

The only endogenous mechanism for doublet crater production that seems 
possible is the formation of calderas. However, terrestrial calderas are not as 
large as most of the Martian craters, and terrestrial volcano-tectonic depres- 
sions, representative of the largest calderas , are markedly noncircular. In the 
past, it usually was assumed that impact of cosmic debris on any planetary body 
should produce random crater patterns. However, recent studies of probable 
impact craters on Earth contradict this assumption. The Clearwater Lakes 
craters (Figure 1) in Canada, consist of one 32-km diameter crater which is 
almost perfectly tangential to a 22.5-km crater. Both craters are generally 
accepted as impact craters because basement rocks have features characteristic 
of those produced by high-pressure shock waves [ Pence , 1965]. There appears 
to be little doubt that two meteoroids impacted simultaneously to produce the 
Clearwater Lakes craters . Only the mechanism leading to the presence of two 
separate impacting bodies above the impact site at the same time is in question. 

Sekiguchi [1970] has recently proposed a theory, which as qualified and 
supplemented in this study, accounts for the splitting of an impacting meteoroid 
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as it approaches a planet or satellite. He concludes that gravitational forces 
acting on an impacting meteoroid may cause internal stresses that exceed the 
meteoroid's tensile strength, causing it to split into bodies of approximately 
equal size. This type of calculation can be compared to the well-known calcula- 
tion of the Roche limit, the critical altitude for which the difference in gravita- 
tional attraction between the planet and a very large body approaching the planet 
is greater than the mutual attraction of the mass elements of the impacting body. 

The Roche limit calculation assumes that the bodies are cohesionless. It is clear 
that bodies producing the Mars craters, with diameters between meters and 
kilometers must have had a finite tensile strength. 

The plausibility of tidal splitting and separation as a means of Mars doublet 
formation will be examined in terms of the required relationships between mass, 
tensile strength, entry velocity, and entry angle of the meteoroid, and the crater 
separation distances and diameters. 

OBSERVATIONS 

Mariner 6 and 7 photographs of cratered terrain in Meridiani Sinus, Deucaleonus 
Regio, and Hellespontus (Figure 2) were studied. Nonrandom effects were first 
noticed for four crater pairs, labeled 1 through 4 in Figure 3. The ordering of 
crater arrangement exhibited by these crater pairs was particularly striking 
because each crater is nearly perfectly tangential to, and of about the same 
apparent age as its neighbor, and craters of a given doublet are nearly equal 


Fig. 2 

Fig. 3 


in size. It seems highly unlikely that doublets with these characteristics could 
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have been formed by random impact of single bodies. Other isolated crater 
pairs are labeled 5 through 8 in Figure 3; crater doublets occurring in groups 
of more than two craters are shown in crater clusters labeled 9 through 13. If 
a doublet is defined as any two craters that overlap or are separated by no more 
than 12.5 km, it can be shown that there are 461 crater doublets present in 
photographs 6N17, -19, -21, -23, and 7N25; this is a very large number when it 
is considered that there are only 906 craters larger than 4 km in the area, 
photographed. 

Figures 4a and b show the size distribution of these craters and the distri- 
bution of overlap or separation for the 461 doublets . The ratio of the radius of 
the largest crater (Rm av) to the radius of the smallest crater (Rmin) of each 
Martian doublet has been determined and the number of Martian doublets as a 
function of Rmax^^min shown in Figure 4c. Inspection of the Figure 4£ and 
the photographs suggests an ordering of crater size for members of doublets 
that is unlikely if all the craters were produced by bodies having equal probability 
of impacting at any position on the surface of Mars . Members of doublets are 
often of similar size. The distribution of Figure 4£ is later compared to similar 
distributions resulting from the Monte Carlo simulations. 

MONTE CARLO SIMULATION MODEL OF MARTIAN CRATERING 

The Monte Carlo model of Martian cratering involves the drawing, at random , of 
craters from the crater size frequency distribution of Figure 4a. Next, random co- 
ordinates are selected for each of the craters drawn from the distribution and these 
are used to determine whether a given crater is close enough to any other crater 
to form a doublet. Random numbers between 0 and 1 are produced by one of two 
pseudorandom number generators of a class known and described as multiplicative 
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congruent generators [ Fassone , et al ., 1968]. That is, the generators are defined 
by the recurrence relationships, Uj+i = aUi- mod(231) where the parameter a is 
defined as 5^ [Marsaglia , 1968] and (2^® + 3) ( IBM 360 Scientific Subroutine 
Package) for the respective generators. 

The selection of a given crater diameter from the size frequency distribution 
of Martian craters (Figure 4a) is accomplished by first drawing a uniform random 
number u. Then, inverse interpolation is performed until the crater diameter D 
is obtained at which the cumulative distribution of diameter lengths equal u. 

Crater diameters given in Figure 4a are actually slightly larger than true values 
because the photographs show slightly oblique views of the surface. However, 
the effect of overestimates of crater diameter is to produce maximum possible 
numbers of doublets in the cratering simulation. An independent calculation 
shows that this leads to overestimates of the number of doublets of less than 
10 percent. Selection of random x,y coordinates for the crater is accomplished 
by generating two numbers between 0 and 1 and multiplying each by the square 
root of the area being cratered (1843.1 km). The roles of the two generators used 
for selection of diameters and coordinates are also reversed in the cratering 
model for checks on randomness. 

The next step in the Monte Carlo simulation is to determine whether a given 
crater, drawn at random, and assigned random coordinates, forms a doublet with 
any of the previous craters. This occurs when any of the following inequalities 
is satisfied; 


+ 
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where is the radius of theith crater drawn and rj is the radius of the jth 
crater, and Xj are the random x coordinates of the Uh crater and jth crater, 
and and_yj are the random y coordinates of theith crater and the jth crater; 

A is the overlap or separation allowed for the particular crater pair in question. 
This value is selected at random from the distribution of A for doublets observed 
on Mars in the same way as crater radius is selected. Each Monte Carlo crater- 
ing simulation involves j - separate solutions of inequality (1) for each of the 
craters 2 . All doublets resulting from satisfactory solutions of these inequalities 
were totaled and recorded by crater number, radius, and coordinates of each mem- 
ber of the doublet. For all doublets, the maximum radius was divided by the mini- 
mum radius and doublets were grouped according to classes of K /R . . The 

° ^ ° —max —mm 

use of all 906 craters in the computation completes one simulation of Martian 
cratering. 

Simulation of the 906 cratering events was repeated 150 times; the resulting 
variation in numbers of doublets was as expected if the Martian cratering were 
random and accomplished 150 times. The statistical characteristics of this dis- 
tribution and the number of doublets by R™„VR^.-„ classes were computed. 

They were then compared to the corresponding Martian data to determine whether 
the cratering is nonrandom and to locate the classes of Martian doublets that con- 
tain most of the doublets produced by the nonrandom process. 

RESULTS OF MONTE CARLO CRATERING MODEL 
The results of the Monte Carlo simulation are shown in Figures 5a and b. 
Figure 5a shows the distribution of the number of doublets that resulted. The 
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average number of doublets that occurred was 137 and two standard deviations 
were 26.2, Thus, there is only 0.3 percent probability that the number of doublets 
should be less than 111 or greater than 163 if craters of doublets arise only from 
random single-body impact. The fact that 461 doublets are observed on the 
Martian surface suggests that most craters are not produced by random single- 
body impact. 

The Monte Carlo result is further substantiated by a separate, analytic 
derivation showing the expected number of doublets produced by random single 
body impact as 141, this quantity is derived as follows: a pair of craters of 
radius rj^, ^2 whose centers are specified on some reference surface by the 
cartesian coordinates (Xj^, Yj^), (X 2 , Y 2 ) form a doublet if the relationship: 


^1 


+ r^ 


A < 


(X, - X^)^ 


+ (Yi 




1/2 


rj^ + rg + A 


(2) 


is satisfied where A denotes the allowable overlap or separation of crater rims. 
This relationship is schematically shown in Figure 6. The area denoted by Aq 
in the figure defines the region within which the center of crater 2 must lie in 
order to form a doublet with crater 1 and is given by: 


Aq = + £.2 + - 7r(Zi + £2 " = 47tA (Ti + X 2 ) (3) 


It is clear that for given values of I 2 ’ probability of two randomly 

distributed craters forming a doublet is given by 


. Ao 4,rA(£^ + r2) 

P (r^, rg, A) = ^ = 


Ar 


Fig. 6 
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where Aj. is given as 3,397 x 10® km. From a probabilistic standpoint, each 

crater radius is between 2.0 km and 125.25 km characteristic of the 906 craters 

)< 

used in this study. On the average, however, two randomly chosen craters are 
expected to have radii equal to the mean radius, denoted here by E(R), of the 906 
craters. Similarly, A may assume any value between 0 km and 12.5 km charac- 
teristic of the 461 crater pairs whose empiric separation or overlap is less than 
or equal to 12.5 km. Thus, again on the average, it is expected that a randomly 
chosen A is expected to have a value equal to the mean A, denoted by AM , of 
the 461 crater pairs. Substituting E(R) for r^^ and and AM for A in equation 
(4) then gives: 


Po = 


_ Sir AM E(R) 


A, 


( 5 ) 


the probability that an arbitrarily chosen pair of randomly distributed craters 
form a doublet. Moreover, since there are possible crater pairs, the ex- 
pected number of doublets among 906 randomly distributed craters is given by 


or 


E(D) = (9®®) 


( 6 ) 


E(D) = 

^ A„ 


( 7 ) 


Estimates of E(R) and AM, respectively, may be obtained from values of the 


arithmetic mean radius of 10.25 km exhibited by the 906 craters used in this 
study and the mean A of 4.53 km exhibited by the 461 crater pairs. Sub- 





of a large central crater. An example is crater 13, which has many satellite 
craters (Figure 2). 


ORIGIN OF NONRANDOM DOUBLET CRATERS AND CRATER CLUSTERS 

Martian Calderas 



resulting from such volcanic activity usually are not circular because they are 
products of collapse and thus reflect local rock inhomogeneities. One property 
of calderas is that the rims are often scalloped due to collapse of lateral pockets 
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features has been observed for the Martian craters. Furthermore, and most 
important, small calderas are not usually circular, and terrestrial calderas equal 
to the median size of the Mars craters are usually volcano-tectonic depressions; 
many do not even resemble Mars craters, for their shapes are controlled by pre- 
existing basement structures and they are often irregular oblong or even rectan- 
gular depressions [ Williams , 1968]. 

From various circularity indexes of terrestrial impact craters and calderas 
and lunar craters Murray [1969] concluded that many lunar craters are of impact 
origin. One circularity index used by Murray is the ratio of largest diameter to 
smallest diameter measured through the crater center of gravity. This index 
was obtained here for Mars craters; isolated craters were measured as well as 
separated members of doublets and clusters. Overlapping craters were not 
included because they would not be circular even if they were produced by the 
impact process. Results of these measurements are shown in Figure 7. Data 
from terrestrial calderas (closed circles) were taken from Murray [1969]; data 
for Martian craters are shown as solid circles. Even with distortion present on 
Martian photographs , which results in overestirhates of the circularity index, the 
mean values for Martian craters are considerably less than those for terrestrial 
calderas. The difference is greatest for large craters because for terrestrial 
calderas departure from circularity increases with crater size. These results 
suggest that the Martian doublets and crater clusters are not calderas— unless 
rock properties are much different on Mars than on Earth, which seems unlikely. 


Fig. 7 
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Crater Doublets Produced by Tidal Splitting of Meteoroids 
The plausibility of tidal splitting and separation as a means of Mars doublet 
formation was examined in terms of the meteoroid mass, tensile strength, entry 
velocity, and entry angle required to produce a crater doublet with crater diam- 
eters Dji^ and D 2 . 

Sekiguchi [1970] shows that tidal forces may split a large meteoroid into two 
fragments of approximately equal size. The fragmentation altitude h in terms of 
the meteoroid radius a and tensile strength T is shown for impact on Mars in 
Figure 8 and is given by: 



( 8 ) 


where p = 3.95 gm/cc is the mean density of Mars and Rq = 3393.4 km is the 
planet radius; a is the meteoroid radius and ^ is the gravitational constant. 
Inspection of Figure 8 shows that if meteoroids have tensile strength characteris- 
tic of hand samples of silicate rocks ( s 10 ® dynes/ cm^) then most of the bodies of 
the size required for forming the Martian doublets (< 20 km diameter) will not be 
split by tidal forces. However, it is probable that large meteoroids are charac- 
terized by tensile strength much less than 10® dynes/cm^. Even hand samples of 
some nieteoroids such as carbonaceous chondrites are very friable and would 
possess tensile strength of less than about lO"^ dynes/ cm^ , However, the strength 
of large meteoroids should be significantly reduced from this value due to the 
presence of large defects or cracks. These would certainly be present if 
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meteoroids have resulted from breakup of a parent body caused by collision 
[ Mason , 1970], Even if the large meteoroids formed as separate bodies, large 
cracks should be present. Differential thermal contraction of different rock forming 
minerals produces stresses which commonly fractures new rock formation on 
earth. 

According to Griffith [1921], there is a basic relationship between compres- 
sive or tensile strength and the length of cracks in any brittle material. The 
crack length can be of atomic dimensions or longer. The theory has been used 
to explain the observed reduction of material strength with increasing size for 
samples small enough to be measured in the laboratory [ Brace , 1961]; it appears 
to be valid for large-scale applications as well. Brace has shown that a modifica- 
tion of the Griffith theory is equivalent to the Coulomb law of failure, which is the 
basis of the Anderson theory of geologic faulting. According to Brace [1961], the 
applied tensile stress at fracture can be calculated from: 

^ /2ES\1/2 

- ~ ( Tra ) 

where E is Young's modulus (dynes/cm^), S is specific surface energy of the 
material in ergs/cm^ and a is the half crack length of the Griffith crack in cm. 

For simplicity, we compute the dependence of tensile strength of rocks as a 
fimction of the hypothetical crack length from equation (9). Following Brace [1961], 
we assume a value of 10^ ergs/cm^ for S of most rocks and a value of 10^^ 
dynes/cm^ for Young's modulus. Figirre 9 shows the calculated relationships 


Fig. 9 


that result. Note that the calculated tensile strength for crack length equal to the 
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grain size of most rocks (without large-scale cracks) is in basic agreement with 
typical observed values of ^ = 10^ dynes/cm2 for hand samples. Meteoroids with 
this tensile strength of the size required to form doublets 20 km) would never 
be split by the gravitational field of Mars. However, any cracks present of the 
order of 10 meters in length (l/lOO or less the size of the body) can reduce the 
tensile strength of the meteoroid to 10^ to 106 dynes/cm^. A meteoroid ^ 20~km- 
diameter with this tensile strength will fail at high altitudes above Mars (Figure 8). 
Thus, Sekiguchi's model for meteoroid breakup [eq. (8)] predicts reasonable values 
of breakup altitude as a function of meteoroid radius and tensile strength for 
meteoroids impacting Mars— providing large cracks are present in meteoroids. 

It may be that if cracks are large enough, T approaches 0 and the altitude of break- 
up would be more accurately approximated by the Roche limit. It is noted that if 
equation (8) is solved for T = 0, we do not obtain the Roche limit. Consideration 
of the relationship between the Roche limit an d Sekiguchi's breakup altitude for 
nearly collisionless bodies is beyond the scope of this paper. However, Sekiguchi's 
result is based on Kelvin's study of tidal forces on the earth in an earth, moon, 
sun system in equilibrium. Strictly speaking a two body system should have been 
developed for the fission model. This model would yield higher breakup altitudes and 
greater numbers of doublets than that reported here. 

Certain conditions of impact velocity and impact angle are required for separa- 
tion of the fission products and formation of craters of the correct diameters for 
production of nearly tangential craters. From the equations of motion of the fission 
products , Sekiguchi shows that the distance of separation of the impact points , Ax, 
on the planetary surface is given by: 
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Ax = Az^cosh ('Y_^ tan 0 (10) 



0 = impact angle from Mars surface normal 
M = 6,455 X 10^® gm = mass of Mars 
R = 3,393.4 km = radius of Mars 

F = 6.664 X 10~® dyne cm^ gm"^ = gravitational constant 
Az = parent impacting body radius 
_t = time of flight (sec) from breakup to impact 
This equation is not the same as that given by Sekiguchi but it is correct. 
Sekiguchi's eq. 5.3 results from a typographical error in the solution for the 
differential equation of motion of the fission products. Equation (10) can be 
modified to yield Ax in terms of crater diameters Di and D 2 produced by the 
fission products of masses Mj^ and M 2 . For impact on this scale, we may assume 
that craters produced by shallow depth of burst explosions simulate impact crater 
formation [ Baldwin , 1963; Oberbeck, 1971]. One of Baldwin's [1963] equations for 
shallow depth of burst explosion craters can be modified to express crater forma- 
tion energy in ergs and crater diameter in cm: 


log ]E = 


log Pj. + 2.924 
0.3284 


( 11 ) 


where E is crater formation energy in ergs and Dg is terrestrial crater diameter 
in cm. Expressing Dg in terms of equivalent diameter on Mars (Dj^) for the same 
impact energy, we obtain from Johnson's [1970] experimental results: 
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5.E = -M( i f=\ where gj^ and are the acceleration of gravity on Mars and 


Earth, respectively. Substituting this value into equation (11) and taking antilogs 
yields in cgs units: 


E = 5.59 X 108 d|^045 


( 12 ) 


The relationship is valid only for vertical impact where all available energy is 
used for cratering. There is evidence that the energy of formation of cratering 
is given by the product of the cosine of the impact angle measured from the sur- 
face normal and the projectile kinetic energy [ Moore , 1971]; we modify equation 
(12) to take account of impact angle and equate energy of formation to the product 
of cosine of impact angle and projectile kinetic energy to yield meteoroid mass M: 

M = 11.18 X 108 p3-045/(v2 cos 6) ( 13 ) 


where V is impact velocity and 0 is impact angle from the Mars surface normal. 
This expression and the assumption of spheroidal shape and p of meteoroid = 3.4 
gm/ cc for the parent body can be used to yield a relationship between parent body 
radius and doublet crater diameters and P 2 that would be produced from 
the fission products of breakup of a parent body of radius Az: 


Az = 


4.32 X 10^ 


V 


2/3 ^osl/3 


e 


( 14 ) 


Substitution for A^ in equation (10) yields Ax, the separation distance of centers 
of the craters in terms of P 2 , V, 0 , M, R and flight time U 
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_ 4.32 X 10^ cosh (y_^ 

cos^/3 0 cot 6 


Finally, division of both sides of this equation by (D^ + P2)/2, the average crater 
diameter, yields a separation ratio y, the ratio of distance between crater centers 
to average crater diameter: 


8.64 X 10^ cosh (yt) 


’ cos 6 cot 6 (Pj + P 2 ) 


Equation (16) can be used to determine the relationship between impact angle 
of the parent body and ratio of separation between crater centers to average crater 
diameter for craters produced at the same impact velocity from the fission prod- 
ucts of breakup at a given altitude. This is accomplished by first measuring D| 
and P 2 of craters forming a given doublet. These values are substituted into 
equation (16), and 9 is varied for constant values of V to compute the separation 
ratio as a fimction of impact angle 6 . Flight time _t , required in equation ( 16 ) 
for each solution for y, is computed from each value of V, 0 , and altitude of 
meteoroid breakup. Altitude of breakup is computed from equation (8) , which 


radius a is equal to Az where Az is given by equation (14) and the values 

assumed for V and 0 ; T is calculated from equation (9) by assuming that the 

largest crack length is always .02 a. Figure 10a shows these solutions Fig. 10 


for the largest isolated observed crater doublet where = 68,5 km and 
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occurs primarily for high entry angles, 9 . Similar solutions have also been obtained 
for doublets with smaller average diameters and for different assumed ratios of 
defect length to meteoroid diameter. The fact that doublet formation is restricted 
to bodies impacting at high angles to the surface normal seems at first consider- 
ation incompatible with the observation that of the total of 906 craters there are 
461 doublets. However, many doublets occur in the same crater cluster and could 
have resulted from breakup of a single body. To more accurately approximate the 
number of impacting bodies that could have been broken up, doublets and clusters 
have been counted as single breakup events. If isolated craters are also counted 
as single impact events, we find that as small a fraction as 3/10 of the impacting 
bodies could have broken before impact. This result is in better agreement with 
the theoretical results of Figure 10, but it is still higher than expected. Some of 
the doublets or clusters counted as single events could have been produced by a 
number of separate impacts. Also, the theoretical model does not consider the 
effect of stress energy stored in the meteoroid or meteoroid rotation. Both effects 
would yield greater separation between fission products and more doublets. Thus, 
the observed number of doublets and clusters seems to be in fair agreement with 
the theoretical predictions of the fission model. 

Estimates of the possible aerodynamic lifting forces were made and included 
in the calculation of tidal separation. In obtaining an estimate of the maximum 
separation due to aerodynamic lift, we ignore drag forces and make the simplify- 
ing assumptions that the velocity v and trajectory angle 6 as measured to the 
normal are constant. According to Berman [1961], the lifting force normal to the 


trajectory is: 
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F2 = 1/2 C_gSp (17) 

where p is the atmospheric density, is the lift coefficient, and S is the 
frontal cross section of the body. Then, we make the approximation 
h = h_o - V cos ® 1 or ^ = -^/v cos 6 , where h denotes the altitude, _t is the 
flight time, and ^ denotes the altitude at _t = 0. Moreover, approximating the 
relationship between p and h [ Berman , 1961] by 

P = PQ±^~ ( 18 ) 


where ^ = 9.86 x 10 and Pq = 1.33 x 10 ^ gm/cc [ Levin et al., 1968], and 
performing a double integration, we obtain the displacement on the surface for a 
body as 


-/3ko 


^0 


= e 


Im cos^0/3^^ 


/3(ko 


“b 


^ - 1 ) 


( 19 ) 


where M = meteoroid mass, S = meteoroid cross section, and i^ = altitude of 
breakup. We then estimate the maximum separation between two meteoroids as 
2 Xq. Separation due to tidal splitting and aerodynamic lift is shown as a solid line 
for each velocity in Figure 10. Aerodynamic lift has little effect on separation due 
to the tenuous Martian atmosphere. 

DISCUSSION 

It is highly probable that a nonrandom cratering effect has contributed to the 
production of the cratered terrain on the planet Mars. It is conceivable that some 
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of the isolated nonrandomly produced doublets could have been produced by 
simultaneous impact of gravitationally coupled asteroids. Tanner [1963] has 
shown that such coupled bodies could remain stable as they impact the Earth. 

Since it has been shown in this study that meteoroids approaching Mars could 
produce the observed doublets, if meteoroids are large enough, meteoroid fission 
due to stresses induced by the gravitational field of Mars is favored at this time 
for the nonrandom production of doublets or crater clusters. 

K the Mars doublets have been produced by simultaneous impact of meteoroids, 
it is reasonable to ejq)ect the development of unusual features due to close prox- 
imity of formation of one crater relative to the other. Doublet craters have been 
produced in sand targets in the laboratory both by separate and simultaneous 
impact of two halves of a cylinder that has been partially severed before launch. 

Complete severing and separation of the cylinder halves occurs as a result of 
spinning induced by rifling in the gim barrel. For doublets produced by simul- 
taneous impact, both pieces of the cylinder are allowed to impact. For doublets, 
produced by separate impact , one-half of the cylinder is captured and one crater 
is produced. Then, on the next firing the other half is caught and the other doublet 
member is produced. Degree of overlap is specified by the ratio of the separation 
between crater centers to diameter of either member crater (S/D). Figure 11 Fig. 11 

shows six doublets produced experimentally; three were formed at the specified values 
of s/d by simultaneous impact of half cylinders and three were formed by separate impact. 

All doublets were produced at about the same impact velocity of 2 .3 km/ sec. One 
interesting result is that the crater wall of one member of a doublet is breached 
by the other member only for those doublets whose members have been formed 
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separately. This is illustrated by the doublets formed simultaneously and 
separately at nearly equal values of S/D = 0.81 and 0.72, respectively. The last 
crater to form the doublet produced by separate impact partially destroyed and 
distorted the rim of the first crater. However, the crater rims of members of 
the doublet produced by simultaneous impact are equally developed. Another 
interesting feature of doublets produced by simultaneous impact (S/p = 0.81) is 
that the common crater wall is straight and very high. A continuation of this wall 
projects onto the target surface so that there is a ridge perpendicular to the 
bilateral axis of symmetry of the doublet. Observations of growing craters shows 
that this is the product of focusing of the ejecta and extreme target deformation. 

For the doublet where the members are widely separated (S/D = 1.3), the ridge 
occurs as a delicate subdued structure exactly midway between the separated 
members. Although exact modeling of the very large Mars doublets is not as- 
sumed or expected, the observations of resulting structures applied to Mars 
craters may support the theory of impact origin. Most isolated Mars doublets 
are characterized by craters that are equally developed, and the common walls 
of many doublets are straight and raised, as in doublets 4 and 5 in Figure 2. This 
may be of special significance because rims of isolated craters usually are not 
well developed on Mars, suggesting strong erosive processes. Focused secondary 
crater fields and ridges are present perpendicular to the bilateral axis of symmetry 
of the doublets produced in the laboratory, but these have not been observed for 
Mars doublets. However, focused secondary crater fields and ridges may have 
been subdued by erosion. It may be possible to observe remnants of such ridges 
and focused ejecta fields on higher resolution photographs that will soon be available. 
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The results of this study are of considerable importance to current and re- 
cent interpretations of the geologic history of Mars and the Moon based on obser- 
vation of Mars cratering density compared to that of the Moon. It has been 
assumed [ Murray , 1971] that Mars craters are the results of random single-body 
impact, contrary to the results presented here. On the basis of this assumption 
and the observed smaller ratio of large Martian craters to small craters on Mars 
as compared to the Moon, Murray concludes that the smaller bodies impacting the 
Moon could be from a different population than those impacting Mars and that the 
paucity of large fresh craters could be the result of an unusually efficient erosive 
process. Both hypotheses may not be true. The observations might be explained 
by the fact that fission of meteoroids appears to be common for impact on Mars, 
and this effect would certainly be less pronounced for impact on the Moon because 
of its smaller gravitational field. Thus, the observed deficiency of large fresh 
Mars craters compared to the Moon could be easily explained by the fact that 
more large bodies may be broken up to produce small craters on Mars than are 
broken up for impact on the Moon. Interpretations based on crater frequency 
alone are only valid when the crater distributions are random and when impact 
mechanisms are the same. 

If the occurrence of nonrandom crater doublets on the surfce of Mars is ex- 
plainable by the fission model, a new and unexpected aspect of the mechanics of 
planetary impact cratering has been revealed by the Mariner photographs. 
Meteoroids larger than some limiting diameter, dependent on the amplitude of 
the gravitational field of the planet and meteoroid tensile strength, might be 
broken up before impact. Thus, conditions on one planet may be more favorable 
to doublet and crater cluster development than on other planets. Other factors such 
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as the role of aerodynamic effects may also influence the number of doublets 
produced. 

There are at least three possible terrestrial doublets that may have been 
produced by splitting due to stresses induced in impacting meteoroids by the 
gravitational field of the Earth. The Clearwater Lakes craters have already been 
described. The Rieskessel and the Steinheim Basin may be members of another 
doublet that are more widely separated than the members of the Clearwater Lakes 
craters. The Odessa and Meteor craters may constitute another widely separated 
pair. 

Finally, we can tentatively suggest that the puzzling concentration of large 
mare impact basins on the Earth-facing side of the Moon, a pronounced large- 
scale example of a nonrandom effect, might in fact be the product of breakup of 
a large impacting body. Such a large body would certainly contain defects large 
enough to reduce the tensile strength to a point where the body is essentially 
cohesionless. In that case, breakup would occur near the Roche limit (2.44 times 
the Moon's radius). Alternatively, the body may have passed close enough to the 
Earth to have been ruptured by the earth's gravitational field. The products of 
breakup might well be a cluster of large impact basins with slight age differences 
that appear greater due to the domination in growth of one basin over others. 
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FIGURE CAPTIONS 

Photograph Clearwater Lakes Craters. Photograph courtesy of 
M. Dence. 

Mars chart showing location of photographs. 

Mars photographs with doublets outlined and labeled. 

Cumulative size frequency distribution of craters on photographs 6N17, 
-19, -21, -23, and 7N25. 

Cumulative number of doublets with overlap or separation distance 
greater than A , km. 

Histogram of the number of observed doublets as a function of ratio of 
largest member diameter to smallest member diameter. 

Distribution of number of doublets that form for 150 Monte Carlo 
simulations of Martian cratering. 

Average and standard deviation of number of doublets formed in Monte 
Carlo simulation with indicated ratio of (solid lines). 

Same distribution for observed doublets, dashed line (from Figure 3c). 
Region within which the center of the second crater must be in order to 
form a doublet in the Monte Carlo simulation. 

Circularly indices of terrestrial calderas (closed circles) from 
Murray [1969], and circularity indices of Martian craters , 

Relationship between breakup altitude above the planet Mars and 
meteoroid radius for different assumed values of meteoroid tensile 


strength. 




- 28 - 


Fig. 9. Calculated relationship between tensile strength and "half crack 
length,” for average silicate rock. 

Fig. 10. Relationship between separation ratio and impact angle for values of 

impact velocity V = 5.15, 10, 20, and 25 km/ sec; (a) tidal effects only 
and (b ) tidal and aerodynamic effects. 

Fig. 11. Crater doublets produced by separate and simultaneous production of 
doublet members. 
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